function belem=localfluidvector(celem)
n=length(celem(1,:));
belem=zeros(2*n-4,1);